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The Russian wheat aphid (RWA), Diuraphis noxia Kurdjumov, is a major global pest of wheat and barley 
production that causes enormous economic damage. Few studies have been conducted to explore and 
decipher the molecular basis of RWA strategies to evade plant defense mechanisms. Gene expression studies 
of RWA in response to wheat genotypes carrying different RWA resistance genes have been initiated in our 
group; however, a secure and accurate understanding of RWA gene expression is dependent on 
identification of suitable reference genes. This study analyzed expression profiles of five potential reference 
genes selected and sequenced during RNA sequencing experiments. The expression of genes coding for actin 
and ribosomal protein L27 was comparatively less variable in RWA fed on different wheat hosts. Results of 
geNorm, NormFinder, and BestKeeper expression analyses support the use of actin and ribosomal protein 
L27 in RT-qPCR studies of RWA gene expression in studies involving RWA-wheat interactions. 



The advent of next-generation sequencing technologies has resulted in a significant increase in transcriptomic 
data for various organisms 1 . The generated transcriptomes have helped researchers not only to decipher 
expression pattern of genes and transcripts but also define the genetic architecture of many species 2 ' 3 . 
Validation of gene expression from such transcriptomic resources has become mandatory for reporting and 
reconfirming expression profiles 4 , and reverse transcription quantitative PCR (RT-qPCR) is rapidly replacing 
traditional methods such as Northern blotting and Ribonuclease Protection Assays (RPA) 5,6 . RT-qPCR's speed, 
sensitivity, efficiency and reproducibility has made it the gold standard for rapid and accurate quantification of 
gene expression profiles 7 from various next- generation sequencing (NGS) datasets. Commercially available 
instruments and consumables have further led to universal acceptance of RT-qPCR 6 . 

Gene expression using RT-qPCR assay is based on the principle of quantifying target mRNA during the 
exponential phase of the PCR. During this phase, the target is doubled with each PCR cycle and the probability 
of PCR-bias due to limited reagents is nullified or decreased. In RT-qPCR, the amplification of product is detected 
on accumulation of fluorescent signal. The cycle at which the fluorescent signal exceeds background signals is 
referred to as the threshold cycle, or Ct (also referred to as the quantification cycle, or Cq). Analysis of Cq is used 
to estimate expression of the respective genes. Several factors, such as RNA quality and quantity, mastermix 
components used in the PCR and efficiency of PCR reaction, influence Cq values 8 . Absolute and relative quan- 
tification methods are generally used to estimate gene expression using RT-qPCR. The absolute quantification 
approach necessitates generating the standard curve of known copy numbers of each target, which in turn 
requires standard curves for multiple targets in a study and knowledge of copy numbers of each target, thus 
limiting its usage. The most widely adopted approach, relative quantification, is based on estimation of gene 
expression normalized to the expression of a control gene known as a reference; therefore, reliable determination 
of reference/references is the central factor in accuracy of this method. 

A gene can be used as a reference when it is highly and stably expressed in all samples under investigation and is 
not co-regulated with a target gene 910 . Reference genes traditionally have been housekeeping genes believed to 
possess inherent stable and constitutive expression irrespective of physiological conditions in different samples or 
treatments under investigation 1112 . The universal stability of housekeeping gene expression was disproven in 
recent years 1314 , negating their use 15 . According to the Minimum Information for Publication of Quantitative 
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Table 1 


Primers designed for amplification and sequencing 


of selected D. noxia genes 
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RPL9 


Acyrthosiphon pisum 39S ribosomal protein 


F:TCGAGGTAATCGATCTGCAA 


55 


245 




L9-LOC1 001 59528 (4e-l 14) 


R : GC ATG AATAAAG ATG AGCCTTG 






RPL5 


Acyrthosiphon pisum 60S ribosomal protein 


F : AGCGTTGGTC AATTTCTTGG 


58 


448 




L5-LOC 100 162441 (0) 


R :TGAAGGTCC AAATGCTTTCC 






TEF2 


Acyrthosiphon pisum elongation factor 


F : C AAGGTATG AATTG ATGTCTGGT 


58 


220 




2-LOC100168198 (2e-86) 


R:TCGTTAAAGCATACTTGCCTGT 







Real-Time PCR Experiments (MIQE) guidelines 16 , reference gene(s) 
are now selected based on their specificity in interactions between a 
species or cell type/s subjected to different treatments or conditions. 

The Russian wheat aphid (RWA), Dim aphis noxia Kurdjumov, 
originally from Central Asia, was introduced globally in the 1900s 17 
and has become a major destructive pest of wheat and barley and 
caused huge economic losses. Host-plant resistance is the most 
acceptable and ecologically beneficial aphid pest management strat- 
egy in wheat growing regions, yet host-plant resistance is transitory 
in many cases, because virulent populations develop to survive in 
once cultivated resistant varieties 18 . The short-term solution to this 
problem traditionally has been to identify new sources of resistance 
and breed new RWA-resistant wheat varieties, but new technologies 
such as NGS and RT-qPCR offer the possibility of understanding the 
molecular principles supporting RWA virulence and engineering 
resistant plants with greater longevity. Our research is focused on 
investigating the molecular mechanisms of RWA-plant interaction 
with the long-term goal of finding novel and durable solutions to 
RWA management. 

The objective of this investigation was to identify a robust RWA 
reference gene for use in validation of gene expression studies in 
RWA-wheat interactions. We identified four previously unreported 
RWA sequences commonly used as reference controls in other bio- 
logical systems and report the most suitable reference controls for 
RT-qPCR assays of RWA genes expressed in aphid- wheat 
interactions. 

Results and Discussion 

Analysis of RNA quality, sequencing and sequence analysis. 

Bioanalyzer analysis of RNA revealed a ribosomal shift due to 
denaturation of 28 s rRNA (see Supplementary Fig. 1 online), 
consistent with previously published examples of insect RNA 19 . 
Multiple peaks were absent, indicating no RNA degradation. Thus, 
screening total RNA using NanoDrop spectrophotometry, gel 
electrophoresis and capillary electrophoresis produced high-quality 
RNA for RNAseq library preparation. The libraries analyzed using 



the Bioanalyzer revealed fragment sizes of >230 bp. The quality of 
the 18 libraries (6 treatments X 3 replications) is shown in 
Supplementary Fig. 2 online. Libraries were sequenced and filtered 
to remove low-quality sequences, providing >10 million reads for 
each library. Results of BLAST analysis and differential gene/ 
transcript analysis will be reported separately. 

Gene selection and RT-qPCR. The commonly used reference genes 
actin, ribosomal protein L9, ribosomal protein L27, transcription 
elongation factor 2 and ribosomal protein L5 were selected based 
on a literature survey. Actin has been used as a reference gene in 
many studies of model insects such as the fruit fly, Drosophila 
melanogaster Meigen 20 . Transcription elongation factor genes were 
evaluated to check their suitability as reference genes in several 
insects, particularly Hemipterans such as the sweetpotato whitefly, 
Bemisia tabaci (Gennadius) 21 and brown plant hopper, Nilaparvata 
lugens (Stal) 22 . Genes coding for ribosomal proteins have been 
reported as most reliable reference in insects such as N. lugens 22 . 
Many transcripts coded for these genes were identified from 
RNAseq data, revealing >93% homology to the respective genes 
in the pea aphid, Acyrthosiphon pisum (Harris) (Table 1). 
Representations of these genes observed in RNAseq data 
confirmed that there was little or no change in their expression. 
These genes were amplified from cDNAs, yielding amplicon sizes 
ranging from 220 bp to 553 bp (Table 1). These sequences have been 
deposited in DNA Databank of Japan with accession numbers 
AB914563-AB914566. 

Amplicon size and melting temperature of RT-qPCR primers 
designed from these sequences are shown Table 2. Primer efficiency 
ranged from 91.50 to 102.65, with corresponding R 2 values from 
0.9922 to 0.9995) (Table 2). Primers were redesigned and efficiencies 
rechecked in cases where the efficiency or R 2 values differed from the 
range of recommended values. Melt curve analysis of all genes 
showed no primer dimers or nonspecific product amplification. 

Analysis using geNorm software. Cq values derived from RT-qPCR 
assay for expression stability were logarithmically transformed as 



Table 2 


Primers designed for amplification and RT-qPCR of 


selected genes from D. noxia 


feeding on different wheat genotypes 


Primer 


Sequence (5'— 3') 


Tm (°C) 


Amplicon size in cDNA (bp) 


RT-qPCR efficiency (%) 


Correlation coefficient (R 2) 


RTActin 


F :TGCCTGATGGTC AAGTC ATC 


55 


110 


91.99 


0.9992 




R :TC AAGGAATGCTTTGAGCTG 










RTRPL27 


F : ACC AGC ACGATTTTACC AGATTTC 


55 


90 


91.50 


0.9995 




RCGTAGCCTGCCCTCGTGTA 










RTRPL9 


F:CATTCCACTGGCAC 1 1 1 1 ICTTT 


55 


100 


92.60 


0.9993 




R:TGCCAGAACGCCCAATTAGT 










RTRPL5 


F : AACCC AAAATGTGTTTACGATGAG 


55 


110 


102.65 


0.9922 




R : AGCTGTCGATGGAGGTTTC AAC 










RTTEF2 


F :C AC ATGGGTC ACCTGGGAATA 


55 


96 


93.07 


0.9993 




R :TCGGATTC ACTGCTGATTTGC 
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Figure 1 | The average expression stability values (M) derived using geNorm software for the candidate genes actin, ribosomal protein L5 (RPL5), 
ribosomal protein L9 (RPL9), translation elongation factor 2 (TEF2) and ribosomal protein L27 (RPL27). The least stable candidate gene is plotted on 
the left, and the most stable gene is on the right. 



input for geNorm analysis. The principle behind this algorithm is 
that if two genes are stably expressed in a sample set then the ratio of 
their logarithmic transformed expression should be constant; thus, 
geNorm ranks genes based on their average expression stability (M), 
and the candidate gene possessing the lowest M value is the most 
stably expressed gene in that set. Figure 1 depicts the average 
expression stability of actin (M = 0.210) and ribosomal protein 
L27 (M = 0.267), the two most stable genes that should be used as 
references in RT-qPCR assays involving D. noxia -wheat interactions. 

geNorm also predicts an optimal number of reference genes for 
accurate representation of gene expression based on calculation of 
pairwise variation (Vn/Vn + 1) between sequential normalization 
factors (NFn and NFn + 1). A cutoff value of 0.15 is considered for 
the ratio, below which there is no requirement of any other reference 
gene 13 . Large pairwise variation represents a significant effect in gene 
expression due to the addition of another gene or genes, reinforcing 
the need for the second gene to be included to derive reliable 



0.75 n 




V2/3 V3/4 V4/5 



Figure 2 | geNorm pairwise variation (V) analysis to determine optimal 
number of reference genes for normalization in RT-qPCR reaction. 



normalization factors. Candidate genes are added based on the rank- 
ing derived using the M value. This analysis, depicted in Figure 2 for 
our experiments, clearly showed that adding a third gene does not 
increase the ratio by more than 0.15; therefore, actin and ribosomal 
protein L27 represent the required reference genes for accurate 
estimation of D. noxia gene expression in aphid-wheat interactions. 

Analysis by NormFinder software. NormFinder analyzes candidate 
reference genes according to inter- and intra-group variation in 
expression. Similar results were obtained using NormFinder, which 
predicted stability value of 0.735 for ribosomal protein L27 followed 
by stability value of 0.738 for actin (Figure 3). Therefore, both 
geNorm and NormFinder outputs provide proof for stable 
ribosomal protein L27 gene that can be used as a reference in all 
the RT-qPCR assays that studies expression patterns of genes in 

1.20 r 




0.20 • 

0.00 I 1 1 1 1 1 

RPL5 RPL9 TEF ACTIN RPL27 

Figure 3 | NormFinder stability index for actin, ribosomal protein L5 
(RPL5), ribosomal protein L9 (RPL9), translation elongation factor 2 
(TEF 2), and ribosomal protein L27 (RPL27). Genes with the least stable 
expression are on the left, and the most stable genes are on the right. 
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Table 3 Descriptive statistics based on 
BestKeeper Software 


crossing point (CP) data 


and expression stability of five D. 


noxio reference 


genes calculated by 


Lr tactor 


Actin 




Ribosomal protein 




Transcription 
elongation factor 2 


L27 


L9 


15 
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17.96 
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07 op 
Z/ Jo 


00 
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i y.zo 


Mean Ar CP 


17.96 


18.78 


27.65 


24.66 


19.33 


min CP 


17.28 


18.07 


23.06 


17.66 


17.18 


max CP 


18.87 


19.54 


33.30 


32.95 


21.82 


CP ± SD 


0.51 


0.47 


3.72 


6.24 


0.97 


CP CV (%) 


2.85 


2.48 


13.47 


25.29 


5.00 


CP CV (%) ± SD 


2.85 ± 0.51 


2.48 ± 0.47 


13.47 ±3.72 5 


± 1.95 


25.29 ± 6.24 



Abbreviations. CP - Crossing point, Geo - geometric, Ar - arithmetic, SD - standard deviation, CV - percentage coefficient of variation, 
n = 1 8 samples for each gene. 



the Russian wheat aphid feeding on different wheat varieties 
containing aphid resistance genes. 

BestKeeper analysis software. BestKeeper considers the Cq values of 
all candidate reference genes, in order to calculate standard deviation 
(SD) and coefficient of variation (CV). The software excludes genes 
with Cq value SDs greater than 1. Results of geNorm and 
NormFinder analyses were further verified when Cq data were 
analysed using BestKeeper software, which predicts gene stability 
based on low CV and SD. Ribosomal protein L27 (CV = 2.48 ± 
0.47%) and actin (2.85 ± 0.51%) were found to be the most stable 
(Table 3). 

Validation of reference genes. Results of the current study indicate 
actin and ribosomal L27 to be the most appropriate reference genes 
in RT-qPCR assays involving RWA-wheat interactions. To validate 
their use, we assessed expression of tRNA-Leu, which we previously 
reported to be up-regulated in the gut transcriptome of RWA2 fed 
Dn4 plants 23 . Results of this experiment revealed significant (p < 
0.001) over-expression of tRNA-Leu in RWA2/ Dn4 interactions in 
comparison to other interactions (Figure 4), and very stable CV- and 
M values (Mean CV = 0.168, Mean M = 0.48) for the actin 
- ribosomal L27 combination. In contrast, use of the actin - RPL5 
combination (Mean CV = 0.97, Mean M = 4.2) or the RPL27 - RPL5 
combination (Mean CV = 0.952, Mean M = 3.95) - resulted in huge 
changes in expression, and Mean CV- and M values above acceptable 



stability values [homogenous (CV < 0.25; M < 0.5) and 
heterogeneous (CV < 0.5; M < 1)]. Thus, the results of the tRNA- 
Leu expression validation experiment confirm the use of actin and 
RPL27 for RWA expression analysis. 

Conclusions. This study emphasizes the use of an appropriate 
reference gene(s) for RT-qPCR studies, reports sequences of five 
selected genes from the Russian wheat aphid, and identifies 
reference genes with stable expression in RWA feeding on wheat 
plants carrying different RWA resistance genes. Analyses of our 
data using three different statistical techniques indicated that RWA 
ribosomal protein L27 (RPL27) is the most stable reference gene, 
closely followed by actin. Ideally, both RPL27 and actin genes 
should be used as references for normalizing expression profiles in 
RT-qPCR studies. This study is an important step in our ongoing 
RNAseq studies to define differential expression of RWA genes in 
RWA-wheat interactions. Results of these experiments will be of 
utmost importance to RWA genome sequencing efforts as well as 
the huge amount of expression data derived from current studies 
using NGS technologies. 

Methods 

Insect and plant material. RWA1 (bio type 1) collected from wheat fields near Hays, 
KS, and RWA2 (biotype 2) individuals collected from wheat fields near Briggsdale, 
CO (via the USDA-ARS Plant Science Research Laboratory at Stillwater, OK), were 
cultured in separate housing in the greenhouse at Kansas State University on plants of 
susceptible wheat cultivar 'Jagger' for use in all experiments. The identity of each 




tRNA-Leu 



(C) 



30-- 



1 20 
I 



1 10-- 



RPL27+RPL5 



tRNA-Leu 



R1D0 



R1D4 



R1D7 



R2D0 



R2D4 



R2D7 



Figure 4 | Relative normalized expression of the t-RNA-Leu gene using reference gene combinations of: (a) actin and ribosomal protein L27 (RPL27), 
(b) actin and ribosomal protein L5 (RPL5) and (c) RPL27 and RPL5. R1D0 = RWA1 fed wheat plants carrying DnO (no resistance) gene; 
R1D4 = RWA1 fed wheat plants carrying the Dn4 resistance gene; R1D7 = RWA1 fed wheat plants carrying the Dn7 resistance gene; R2D0 = RWA2 fed 
wheat plants carrying DnO (no resistance) gene; R2D4 = RWA2 fed wheat plants carrying the Dn4 resistance gene; R2D7 = RWA2 fed wheat plants 
carrying the Dn7 resistance gene. 
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biotype was verified independently in diagnostic plant differential greenhouse assays 
at Stillwater and Manhattan. Voucher specimen no. 176 (RWA2) is deposited with 
the Kansas State University Museum of Entomological and Prairie Arthropod 
Research. 

RWA1 and RWA2 adults were starved for 12 h, then fed on wheat plants con- 
taining RWA-resistant genes Dn4 (Yumar), Dn7 (94M370), and the susceptible 
control Yuma (DnO) in a fine screen mesh cage in the greenhouse. All plants were 
grown in the greenhouse in 16.5-cm-diameter-plastic pots containing Pro-Mix-Bx 
potting mix (Premier ProMix, Lansing, MI). Environmental conditions were 24:20°C 
day/night and a photoperiod of 14 : 10 [L:D] h). Aphids were collected from plants at 
24, 48, 72 and 96 h post-feeding, stored immediately in RNAlater solution (Qiagen, 
GmbH, Hilden, Germany), incubated overnight at 4°C and frozen to stabilize total 
RNA according to manufacturer's recommendations. 

D. noxia RNA isolation and quality check. Total RNA was isolated from each 
sample (24, 48, 72 and 96 h post infestation) separately. Silica matrix-based RNA 
extraction was performed using the RNeasy Plus Kit (Qiagen). Tissues were 
homogenized using a micropestle according to the manufacturer's recommendations 
for isolation with some modifications. The gDNA eliminator column was not used for 
DNase digestion. Instead, an extra on-column-based digestion was performed using 
an RNase-free DNase set (Qiagen). Purity and integrity of isolated total RNA was 
assessed by three different methods, including absorbance at 260 and 280 nm in a 
NanoDrop spectrophotometer (Thermoscientific, Wilmington, DE), 1% agarose 
(RNase-free grade) gel electrophoresis using GelGreen staining (Biotium Inc., 
Hayward, CA) and capillary electrophoresis using an RNA Nano Lab -Chip (Agilent, 
Santa Clara, CA) and an Agilent 2100 Bioanalyzer system 24 . 

Total RNAs qualifying according to recommended NanoDrop spectrophotometry 
quality standards were subjected to agarose gel electrophoresis, and samples with 
intact bands were processed in the Agilent 2100 Bioanalyzer at the Kansas State 
University Integrated Genomics Facility in Manhattan, KS. The Bioanalyzer provided 
details about RNA integrity summarized as a RNA integrity number (RIN) (from 1-- 
10; maximum-no degradation) and rRNA (28 s/18 s) ratio (~2 means very high 
quality RNA) 25 . Total RNA collected at 24, 48, 72 and 96 h post-feeding by each 
biotype was pooled to form the treatments RWA1 fed DnO (KW Al/ DnO) y RWA1 fed 
Dn4 (RWAl/0n4), RWA1 fed Dn7 (RWA1/Dn7), RWA2 fed DnO (RWAl/DnO), 
RWA2 fed Dn4 (RWA2/Drc4), and RWA2 fed Dn7 (RWA2/Drc7). 

TrueSeq library construction, sequencing and bioinformatics analysis. mRNA was 
purified, using oligo-dT beads, from 1 jag of total RNA of each treatment- replication 
sample provided in TruSeq RNA Sample Preparation Kit v2 (Illumina Inc., San 
Diego, CA). After purification, first and second strand cDNA was synthesized, 3 ' ends 
were repaired and adenylated, and DNA fragments were enriched after adapter 
ligation and indexing. The generated library was validated and pooled based on the 
specific index or bar-code for sequencing. Samples were loaded in two lanes and 
sequenced in a Hiseq 2500 Illumina sequencer at the University of Kansas Medical 
Center in Kansas City, KS. 

Raw sequences generated from sequencing were de-multiplexed to define their 
sample of origin. Fastq output files of each sample were assessed for quality control 
(QC) using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), 
and de novo assembly of sequences was conducted using Trinity software (http:// 
www.trinity-software.com). Transcript and gene abundance were estimated and 
expression was evaluated using DESeq (differential gene expression analysis) and 
EdgeR (empirical analysis of digital gene expression), both from Bioconductor Open 
Source Software for Bioinformatics (http://bioconductor.org/). 

cDNA synthesis and selection of candidate reference genes. cDNA used for PCR 
standardization and gene amplification was synthesized using Superscript III First- 
strand Synthesis Supermix (Life Technologies, NY) and oligodT primers according to 
manufacturer's protocol. For RT-qPCR, an equal amount of total RNA (500 ng) from 
each sample was used for cDNA synthesis, using the iScript Reverse Transcription 
Supermix for RT-qPCR (Bio-rad, Hercules, CA). 

Candidate reference genes were selected based on literature survey and iden- 
tification in the library. Primers were designed using Primer 3 (http://primer3.ut.ee/) 
(Table 1). The PCR profile used was 94°C for 30 s, 30 cycles of 94°C for 30 s followed 
by 55-58°C and 72°C, each for 30 s. The amplicons were sequenced in the Kansas 
State University Integrated Genomics Facility. 

RT-qPCR primer design and test of amplification efficiency. Primers of sequenced 
transcripts considered for RT-qPCR were designed using Primer Express Software 
V.3.0.1 (Life Technologies) with slight modifications of the default parameters 
(Table 2) to avoid hairpin, self-dimer and cross-dimer secondary structures. Primers 
were selected based on low penalty scores, which indicated a good match to the set 
parameter settings. 

Amplification efficiency of primers was estimated by SYBR Green chemistry RT- 
qPCR assay using 1, 10~\ 10" 2 , 10" 3 and 10" 4 -fold dilutions of pooled cDNAs ofthree 
technical replicates for each gene. Each PCR mixture consisted of 1 ul of cDNA, 
250 mM of primer and 12.0 ul of iTaq Universal SYBR Green Supermix (Bio-Rad 
Laboratories) in a total volume of 20 ul. The RT-qPCR profile used was 95°C for 
3 min followed by 40 cycles of 95°C for 10 s, 55°C for 30 s and 72°C for 30 s. After 
RT-qPCR, melt curve analysis was also performed to check for unwanted products or 
contaminants during the PCR reaction. Three technical replicates were used for 
standard curve generation. The fluorescence signals were captured, quantification 



cycle (Cq) values were tabulated and each mean Cq was plotted against the logarithm 
of the pooled cDNA dilution factor. The slope in a linear regression model was used to 
determine the efficiency (E = [io (1/_slope) — 1] X 100%) for each primer. Correlation 
coefficient (R 2 ) values were calculated using the standard curve. 

RT-qPCR for calculation of expression stability and statistical analysis. The 

expression stability of the selected genes was calculated in RT-qPCR assay with 
selected primers and PCR conditions stated above. Each PCR mixture constituted of 
1 ul of cDNA from each treatment sample, 250 mM of each primer, and 12.0 ul of 
iTaq Universal SYBR Green Supermix (Bio-Rad Laboratories) in a total volume 
of 20 ul. The RT-qPCR profile was 95°C for 3 min followed by 40 cycles of 95°C for 
10 s, 55°C for 30 s and 72°C for 30 s. Each reaction was performed for three technical 
replicates and three biological replicates in 96-well optical grade PCR plates sealed 
with optical sealing tape (Bio-Rad Laboratories). Melt curves were generated for each 
sample by heating the PCR amplicon from 65°C to 95°C with a gradual increase of 
0.5°C per cycle of 0.5 s. 

The Microsoft Excel-based software geNorm 13 , NormFinder 26 and BestKeeper 27 
was used to rank reference genes for stability of expression across all experimental 
samples. For geNorm and NormFinder analyses, Cq values were transformed to 
relative quantities using the formula 2" ACq , where ACq = Cq of the gene in selected 
sample - minimum Cq of corresponding gene in the experiment. The sample with 
minimum Cq or maximum expression was used as the calibrator with a set value of 1. 
No transformed Cq values are required for BestKeeper analysis. Other procedures 
were followed according to the instructions mentioned in the software manual. 



1 . Hudson, M. E. Sequencing breakthroughs for genomic ecology and evolutionary 
biology. Mol. Ecol. Resour. 8, 3-17 (2008). 

2. Adams, J. Transcriptome: connecting the genome to gene function. Nature Edu. 1, 
195 (2008). 

3. Stevenson, B. J., Iseli, C, Beutler, B. & Jongeneel, C. V. Use of transcriptome data 
to unravel the fine structure of genes involved in sepsis. /. Infect Dis. 187, S308-14 
(2003). 

4. Paim, R. M. et al. Validation of reference genes for expression analysis in the 
salivary gland and the intestine of Rhodnius prolixus (Hemiptera, Reduviidae) 
under different experimental conditions by quantitative real-time PCR. BMC Res 
Notes. 5, 128. doi: 10.1186/1756-0500-5-128 (2012). 

5. Klie, M. & Debener, T. Identification of superior reference genes for data 
normalisation of expression studies via quantitative PCR in hybrid roses (Rosa 
hybrida). BMC Res. Notes 4, 518 (2011). 

6. Kheirelseid, E. A., Chang, K. H., Newell, J., Kerinand, M. J. & Miller, N. 
Identification of endogenous control genes for normalisation of real-time 
quantitative PCR data in colorectal cancer. BMC Mol. Biol. 11, 12. doi: 10.1186/ 
1471-2199-11-12 (2010). 

7. Arya, M. et al. Basic principles of real-time quantitative PCR. BioTechniques 5, 
209-219 (2005). 

8. Mahoney, D. J. et al. Real-time RT-PCR analysis of housekeeping genes in human 
skeletal muscle following acute exercise. Physiol. Genomics 18, 226-231 (2004). 

9. Thellin, O. et al. Housekeeping genes as internal standards: use and limits. 
/. Biotechnol. 75, 291-295 (1999). 

10. Lee, P. D., Sladek, R., Greenwood, C. M. & Hudson, T. J. Control genes and 
variability: absence of ubiquitous reference transcripts in diverse mammalian 
expression studies. Genome Res. 12, 292-297 (2002). 

11. Bustin, S. A. Developments in real-time PCR research and molecular diagnostics. 
Expert Rev. Mol. Diagn. 10, 713-715 (2010). 

12. Butte, A. J., Dzauand, V. J. & Glueck, S. B. Further defining housekeeping, or 
"maintenance," genes focus on "A compendium of gene expression in normal 
human tissues." Physiol. Genomics 7, 95-96 (2001). 

13. Vandesompele, J. et al. Accurate normalization of real-time quantitative RT-PCR 
data by geometric averaging of multiple internal control genes. Genome Biol. 3, 
research0034. 1-0034. 11 (2002). 

14. Klie, M. & Debener, T. Identification of superior reference genes for data 
normalisation of expression studies via quantitative PCR in hybrid roses (Rosa 
hybrida). BMC Res. Notes 4, 518 (2011). 

15. Gutierrez, L., Mauriat, M., Pelloux, J., Belliniand, C. & Van Wuytswinkel, O. 
Towards a systematic validation of references in Real-Time RT-PCR. Plant Cell 
20, 1734-1735 (2008). 

16. Nolan, T., Hands, R. E. & Bustin, S. A. Quantification of mRNA using real-time 
RTPCR. Nat Protoc. 1, 1559-1582 (2006). 

17. Morrison, W. P. & Peairs, F. B. [Response model concept and economic impact] 
Response Model For An Introduced Pest-The Russian Wheat Aphid [Quisenberry, 
S.S. & Peairs, F. B. (eds.)] [1-11] (Thomas Say Publications in Entomology, 
Lanham, MD, 1998). 

18. Marasas, C. Wheat management practices and adoption of the Russian wheat 
aphid integrated control programme in the study area. Socio-economic impact of 
the Russian wheat aphid integrated control programme. PhD Thesis, University 
of Pretoria, Pretoria, South Africa. (1999). 

19. Winnebeck, E. C, Millar, C. D. & Warman, G. R. Why does insect RNA look 
degraded? /. Insect Sci. 10, 159 (2010). 

20. Ponton, F., Chapuis, M.-P., Pernice, M., Sword, G. A. & Simpson, S. J. Evaluation 
of potential reference genes for reverse transcription-qPCR studies of 



SCIENTIFIC REPORTS | 4:5059 | DOI: 1 0.1 038/srep05059 



5 



physiological responses in Drosophila melanogaster. J. Insect Physiol. 57, 840-850 
(2011). 

21. Li, R. et al. Reference gene selection for qRT-PCR analysis in the sweet potato 
whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae). Plos One 8, e53006. 
doi:10.1371/journal.pone.0053006 (2013). 

22. Yuan, M. et al. Selection and evaluation of potential reference genes for gene 
expression analysis in the brown planthopper, Nilaparvata lugens (Hemiptera: 
Delphacidae) using reverse-transcription quantitative PCR. PLoS One 9, e86503. 
doi: 10. 1 37 l/journal.pone.0086503 (2014) . 

23. Anathakrishnan, R. et al. Comparative gut transcriptome analysis reveals 
differences between virulent and avirulent Russian wheat aphids, Diuraphis 
noxia. Arthropod Plant Interact 8, 79-88 (2014). 

24. Van de Goor, T. A. (2003) The principle and promise of labchip 
technology .Pharma Genomics 3, 16-18 (2003). 

25. Schroeder, A. et al. The RIN: an RNA integrity number for assigning integrity 
values to RNA measurements. BMC Mol. Biol. 7, 3 (2006). 

26. Andersen, C. L., Jensen, J. L. & Orntoft, T. F. Normalization of real-time 
quantitative reverse transcription-PCR data: a model-based variance estimation 
approach to identify genes suited for normalization, applied to bladder and colon 
cancer data sets. Cancer Res. 64, 5245-5250 (2004). 

27. Pfaffl, M. W., Tichopad, A., Prgomet, C. &Neuvians, T. P. Determination of stable 
housekeeping genes, differentially regulated target genes and sample integrity: 
BestKeeper--Excel-based tool using pair- wise correlations. Biotechnol. lett. 26, 
509-515 (2004). 



Acknowledgments 



The authors thank Drs. Kan Yan Zhu and M-S Chen for their critical reviews of the 
manuscript. This research was supported by the Kansas Wheat Commission, the Kansas 
Wheat Alliance and Kansas State University. This is contribution 14-240-J from the Kansas 
Agricultural Experiment Station. 

Author contributions 

D.K.S. and C.M.S. scoped and designed the study; D.K.S. performed the experiments and 
analyzed data; D.KS. and C.M.S. interpreted results and wrote the manuscript. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Sinha, D.K. & Smith, CM. Selection of Reference Genes for 
Expression Analysis in Diuraphis noxia (Hemiptera: Aphididae) Fed on Resistant and 
Susceptible Wheat Plants. Sci. Rep. 4, 5059; DOI:10.1038/srep05059 (2014). 



This work is licensed under a Creative Commons Attribution-NonCommercial- 
NoDerivs 3.0 Unported License. The images in this article are included in the 
article's Creative Commons license, unless indicated otherwise in the image credit; 
if the image is not included under the Creative Commons license, users will need to 
obtain permission from the license holder in order to reproduce the image. To view 
a copy of this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



SCIENTIFICREPORTS | 4:5059 | DOI: 1 0.1 038/srep05059 



6 



